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Abstract 

We present a theoretical study of Weiss oscillations in magnetoconductivity of bilayer graphene. 
Bilayer graphene in the presence of a perpendicular magnetic field and a unidirectional weak elec- 
tric modulation is considered. We determine the a yy component of the magnetoconductivity tensor 
for this system which is shown to exhibit Weiss oscillations. We show that Weiss oscillations in the 
magnetoconductivity of bilayer graphene are enhanced and more robust with temperature com- 
pared to those in conventional two-dimensional electron gas systems whereas they are less robust 
with temperature compared to monolayer graphene. In addition, we also find phase differences of it 
and 2tt in the magnetoconductivity oscillations compared to monolayer graphene and conventional 
2DEG system which arises due to the chiral nature of quasiparticles in graphene. 
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I. INTRODUCTION 



The successful preparation of monolayer graphene has allowed the possibility of study- 
ing the properties of electrons in graphene The nature of quasiparticles called Dirac 
electrons in these two-dimensional systems is very different from those of the conventional 
two-dimensional electron gas (2DEG) realized in semiconductor heterostructures. Graphene 
has a honeycomb lattice of carbon atoms. The quasiparticles in graphene have a band 
structure in which electron and hole bands touch at two points in the Brillouin zone. At 
these Dirac points the quasiparticles obey the massless Dirac equation. In other words, 
they behave as massless Dirac particles leading to a linear dispersion relation = vk ( 
with the characteristic velocity v ~ 10 6 m/s). This difference in the nature of the quasi- 
particles in graphene from conventional 2DEG has given rise to a host of new and unusual 
phenomena such as anamolous quantum Hall effects and a it Berry phase l| 2J. Besides the 
fundamental interest in understanding the electronic properties of graphene there is also 
serious suggestions that it can serve as the building block for nanoelectronic devices 

In addition to the graphene monolayer, there has been a lot of interest in investigating the 
properties of bilayer graphene. The quasiparticles in bilayer graphene exhibit a parabolic 
dispersion relation which implies that they are massive particles. These quasiparticles are 

n n n n q 

also chiral and are described by spinor wavefunctions|2j, |4], [5|, la Recent theoretical work 
on graphene multilayers has also shown the existance of Dirac electrons with a linear energy 
spectrum in monolayer graphene and a parabolic spectrum for Dirac electrons in bilayer 4j . 
Bilayer graphene consists of two monolayers stacked as in natural graphite. This, Bernal 
stacking, yields a unit cell of four atoms with the result that there are four electronic bands. 
In k space, the bilayer has a hexagonal Brillouin zone. Its physical properties are mainly 
determined by the eigenvalues and eigenfunctions at two inequivalent corners of the Brillouin 
zone, K and K', where the it* conduction and n valence bands meet at the Fermi surface. 
Due to the strong interlayer coupling both the conduction and valence bands are split by 
an energy ~ OAeV near the K and K' valleys. Only two of these bands, upper valence and 
lower conduction band, are relevant at low energy and they can be described by the effective 
Hamiltonian given below 

It was found years ago that if conventional 2DEG is subjected to artificially created 
periodic potentials it leads to the appearence of Weiss oscillations in the magnetoresistance. 
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This type of electrical modulation of the 2D system can be carried out by depositingan a rray 



of parallel metallic strips on the surface or through two interfering laser beams 



io| 



Weiss oscillations were found to be the result of commensurability of the electron cyclotron 



diameter at the Fermi energy and the period of the e 
were found to be periodic in the inverse magnetic field 



ectric modulation. These oscillations 



Recently, an investigation of 



electric field modulation effects on transport properties in monolayer graphene has led to the 



prediction of enhanced Weiss oscillations in the magnetoconductivity 



111 ]. In addition, the 



magnetoplasmons spectrum, density of states, bandwidth and thermodynamics properties 
of monolayer graphene in the presence of electrical modulation have been investigated so 
far [l3j|. In this work we are interested in studying the effects of electric modulation on 
magnetoconductivity in bilayer graphene and we compare the results obtained with those 
of monolayer graphene and the conventional 2DEG. 



II. FORMULATION 

We consider symmetric bilayer graphene within the single electron approximation de- 
scribed by the following effective Hamiltonian (H — c—1 here) 



1 / {P x ~lPy 
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H = -— I " \- ^ — yj I (1) 
2m \ ( Px + lPy f 

where ~j) = -iV — e A , with the vector potential expressed in the Landau gauge as 
A = (0, Bx, 0) and the magnetic field is B = (0, 0, Bz) , which is perpendicular to the 
bilayer graphene, m is the effective mass of the electrons in bilayer: m = 0.043m e with m e 
the usual electron mass. The energy eigenvalues and eigenfunction in the presence of the 
magnetic field are 

e(n) = u c y/n(n-l), n>2 (2) 

where uj c — — is the cyclotron frequency. For the low magnetic fields considered in this work, 
the Hamiltonian of Eq.(l) and the Landau level spectrum in Eq.(2) adequately captures 
the low energy electronic properties in bilayer in the presence of a magnetic fieldQ]. The 
eigenfunction can be written as 



where L y is the y— dimension of the bilayer and the normalized harmonic oscillator eigen- 
function are 

®n{x) = -^===exp l ^ « nH n {——), 

with center of the cyclotron orbit x = l 2 k y . We now consider a weak one-dimensional 
periodic electric modulation in the x— direction given by the following Hamiltonian 

H' = VoCOs(Kx), (4) 

where K = 2n/a , a is the period of modulation and Vq is the amplitude of modulation. We 
apply standard perturbation theory to determine the first order correction to the unmodu- 
lated energy eigenvalues in the presence of modulation with the result 

e'(n,x ) = V n cos(Kx ), (5) 

where 

Vq 

v n{u) = — exp(-u/2)(L n («) +L„_ 2 (m)), 

u = K 2 l 2 /2, and L n (u) are Laguerre polynomials. 

From equations (2) and (5), the energy eigenvalues for the system in the presence of 
modulation are 



e(n, xo) = lq c \J n(n — 1) + V n cos(Kxq). (6) 

From equation (6) we observe that the formerly sharp Landau levels are now broadened 
into minibands by the modulation potential. Furthermore, the Landau bandwidth (~| V n |) 
oscillate as a function of n, since L n (u) is an oscillatory function of its index. 

The bandwidth contains an average of Laguerre polynomials with indices n and n — 2. 
To compare, in the electrically modulated monolayer graphene the bandwidth depends on a 
linear combination of Laguerre polynomials with indices n and n — 1 whereas for standard 
electrons in 2DEG there is only a single term that contains Laguerre polynomial with index 
n. We expect that this modulation induced change in the electronic density of states to 
influence the magnetoconductivity of bilayer graphene and this is calculated in the following 
section. 
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III. MAGNETOCONDUCTIVITY WITH PERIODIC ELECTRIC MODULA- 
TION 



To determine the magnetoconductivity in the presence of weak electric modulation we 
apply the Kubo formula in the linear response regime. In the presence of the magnetic 
field, the main contribution to the Weiss oscillations in magnetoconductivity arises from 

'his is the diffusive conductivity 



scattering induced migration of the Larmor circle center^ 
and we shall determine it following the approach in 1(J 



11 



121 ] . In the case of quasielastic 



scattering of the electrons, the diagonal component a yy of the conductivity can be calculated 
by the following expression, 

°w = E^cM 1 ~ /( £ c)]^c)(*4) 2 (7) 

L x , L y , are the dimensions of the layer, (3 = is the inverse temperature with k B the 
Boltzmann constant, f(e) is the Fermi Dirac distribution function and r(e) is the electron 
relaxation time and ( denotes the quantum numbers of the electron eigenstate.The diagonal 
component of the conductivity a yy is due to modulation induced broadening of Landau 
bands and hence it carries the effects of modulation in which we are primarily interested in 
this work. a xx does not contribute as the component of velocity in the x-direction is zero 
here. The collisional contribution due to impurities is not taken into account in this work. 
The summation in Eq.(7) over the quantum numbers ( can be written as 

£ q n=0 n=0 

where A = L x L y is area of the system. The component of velocity required in Eq.(7) can 
be calculated from the following expression 

d 

Substituting the expression for e(n, x ) obtained in Eq.(6) into Eq.(9) yields 

vl = — — — sm{Kx ). (10) 

With the results obtained in Eqs.(8), (9) and (10) we can express the diffusive contribution 
to the conductivity given by Eq.(7) as 

Oyy = A Q (f) (11) 



where 



A = -V 2 e 2 r{3 (12) 

7T 



and the dimensionless conductivity of bilayer graphene <fi is given as 



(13) 



n=0 



where g{e) = exp[/3(e — Bp] and Bp is the Fermi energy 

IV. ASYMPTOTIC EXPRESSIONS 

To get a better understanding of the results of the previous section we will consider 
the asymptotic expression of conductivity where analytic results in terms of elementary 



functions can be obtained Hj. We shall compare the asymptotic results for the dimensionless 
conductivity obtained in this section with the results obtained for the electrically modulated 
conventional 2DEG system. We shall also compare these results with recently obtained 
results for the monolayer graphene that is subjected to only the electric modulation. 

The asymptotic expression of dimensionless conductivity can be obtained by using the 
following asymptotic expression for the Laguerre polynomials 

1 7T 

exp~" /2 L n {u) — > —= cos(2^/nu - -). (14) 



Note that the asymptotic results are valid when many Landau Levels are filled. We now 
take the continuum limit: 

oo 

»-->^,E--> /- (15) 

to express the dimensionless conductivity in Eq.(13) as the following integral 

oo 

, 1 f , g(e) pa 2 » — — 2 , — 7r 
(P = ^J [g(e) + lW \ln COS (vV n ) cos VVrm-j) (16) 
o 

where u = 2n 2 /b and the dimensionless magnetic field b is introduced as b = with 

ear 

Now assuming that the temperature is low such that (3~ l <C £f and replacing e = 
Ef + s(3~ l , we rewrite the above integral as 



V^ ms 'f^ f , idse ' v> cos* (^-l + ^Ls) (17) 



4/3 \ p J J (e s + l) 2 \ b 4 ptu c 



where p = kpa = \/27m e a is the dimensionless Fermi momentum of the electron. To obtain 
an analytic solution we have also replaced e by £p in the above integral except in the sine 
term in the integrand. 

The above expression can be expressed as 



j2/e F bu c 2 / 2tt\ 
cos — 



ds 



cos 



2irp 7T 47T 



4/3 "" \p J J cosh 2 (s/2) \ b 4 jw c 

— oo 

The above integration can be performed by using the following identity 



dx 



cos ax 



an 



with the result 



4vr 2 T, 



■ cos 



B 



2tt 
P 



cosh 2 /3x 2/3 2 sinh(a7r/2/3) 



i - A[ B +2A (i Uos2 



b 8 



(18) 



(19) 



(20) 



where Tb is the characteristic damping temperature of Weiss oscillations in bilayer graphene: 



k R T, 



bp 



Ait 2 ma? 
bp 



and A(x) 



sinh(x) 



_(x—>oo) _ 2xf} _ 



V. COMPARISON WITH ELECTRICALLY MODULATED MONOLAYER 
GRAPHENE 



We will now compare the results obtained in this work with results obtained in ll|] for 
the case of electrically modulated monolayer graphene system. We will first compare the 
energy spectrum in the two cases. The difference in the energy spectrum due to modulation 
effects was obtained in Eq.(6). If we compare this result with the corresponding expression 
for the electrically modulated monolayer graphene case, we find the following differences: 
Firstly, in the monolayer we have an average of two successive Laguerre polynomials with 
indices n and n — 1 whereas here we also have the average of two Laguerre polynomials but 
not successive ones but rather with indices n and n — 2. Secondly, in the monolayer the 
energy eigenvalues are multiplied by the square root of the Landau band index sfn whereas 



in the bilayer we have \Jn(n — 1) factor. Thirdly, the cyclotron frequency in the two systems 
is different since the quasiparticles in monolayer are massless Dirac particles whereas they 
have a finite mass in the bilayer. These differences cause the velocity expression for the 
electrons given by Eq.(10) to be different in the two systems. 



We now compare the expressions for dimensionless conductivity <fi given by Eq. (20) with 



the electrically modulated case (Eq.(22) in 11|). The argument of the cosine terms in the 



expression for bilayer are 2ir/p whereas in monolayer it is n/p which results in the phase 
difference of n in the the dimensionless conductivity in the two systems. This we expect as 
the quasiparticles in graphene (both monolayer and bilayer) are chiral and acquire a Berry's 
phase in the presence of a magnetic field The Berry's phase acquired by Dirac electrons in 
monolayer graphene is tt whereas it is 2ix for particles in bilayer graphenej^, Sj- Therefore we 
observe a difference in phase of tt in the magnetoconductivity oscillations in the two systems. 
The dimensionless magnetoconductivity for both electrically modulated mono- and bi-layer 
graphene as a function of inverse magnetic field is shown in Fig.(l)at temperature T = 6K , 
electron density n e = 2.3 x 10 11 cm -2 and period of modulation a = 350nm.We also observe 
that in the region of high magnetic field SdH oscillations are superimposed on the Weiss 
oscillations. The oscillations are periodic in 1/B and the period depends on electron density 
as \FnZ. 



VI. COMPARISON WITH STANDARD ELECTRON SYSTEM IN 2DEG 

We start by comparing the energy spectrum and velocity expression obtained in Eq.(6) 
and Eq.(10) with similar expressions for the conventional 2DEG where the the quasiparticles 
are standard electrons . For the energy spectrum, we find that the Landau level spectrum 
is significantly different from that of standard electrons in conventional 2DEG. The first 
term uj c ^n(n — 1) in Eq.(6) has to be compared with uj c (ti + 1/2) with uj c = eB/m e for 
standard electrons. Not only the dependence on the Landau level index n is different in the 
two systems but the cyclotron frequency is also not the same due to the difference in mass of 
the quasiparticles. The modulation effects are carried by the second term where the essential 
difference is in the structure of the function V n (u) = ^ exp(— u/2)(L n (u) + L n _ 2 (w)). We 
find that there is a basic difference: In bilayer we have a sum of two Laguerre polynomials 
with indices n and n — 2 whereas only a single Laguerre polynomial appears in the corre- 
sponding term for standard electrons in 2DEG. This difference in the V n (u) function causes 
the velocity expression for the electrons in bilayer given by Eq.(10) to be different from that 
of the standard electrons. To highlight the difference in the dimensionless conductivity in 
the two systems, we compare the asymptotic expression in bilayer Eq.(20) with the corre- 
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sponding expression for 2DEG (Eq. (25) in [ll|). We find that dimensionless conductivity 
in bilayer has an additional prefactor cos 2 (jjf^J which is not present in the corresponding 
expression for 2DEG. In addition, conductivity in bilayer contains the characteristic damp- 
ing temperature Tg which is higher than the corresponding damping temperature in 2DEG 
T p due to the smaller effective mass of the quasiparticles in bilayer. This results in the 
magnetoconductivity oscillations to be more robust with temperature than in 2DEG. To 
see the effects of this difference on the magnetoconductivity we present the dimensionless 
magnetoconductivity for both electrically modulated bilayer graphene and the electrically 
modulated standard 2DEG in Fig. (2), as a function of inverse magnetic field at temperature 
T = 6K , electron density n e = 2.3 x 10 ll cm~ 2 and period of modulation a = 350nm. We 
find that the there is a difference in phase of 2ir between the oscillations in magnetoconduc- 
tivity in the two systems since the quasiparticles in bilayer graphene are chiral. A Berry's 
phase of 2tt is acquired by the quasiparticles in bilayer relative to the standard electrons 
resulting in the appearence of 2tt phase difference in the magnetoconductivity oscillations. 
We also find a peak in magnetoconductivity in 2DEG that is absent in bilayer which is due 
to the absence of contribution from the n = and n = 1 Landau levels as they lie at zero 
energy. 

We also find that the magnetoconductivity oscillations in bilayer graphene are less 
damped by temperature and are more prounced as compared to those in conventional 2DEG 
system whereas they are less pronounced and are more damped with temperature compared 
to those in monolayer graphene. This can be seen in Fig. (3) where dimensionless con- 
ductivity as a function of inverse magnetic field is presented for the three systems. The 
parameters used are: T = 6K , electron density n e = 2.3 x 10 n cm~ 2 and period of modula- 
tion a = 350nm. This can be understood by considering the temperature scale for damping 
of Weiss oscillations in bilayer graphene obtained from Eq.(20) which is characterized by 
Tr given as HbTb = 4 2 &p 2 whereas the characteristic tempererature for 2DEG is given in 



IfJ 



ll| as ksTp = 4 2 ^ a 2 ■ Comparing T B and T p the essential difference is the difference 



in the effective masses of the quasiparticles in the two systems. Since the quasiparticles 
in bilayer have a smaller effective mass m = 0.043m e , the characteristic damping tempera- 
ture Tg is higher in bilayer than in conventional 2DEG characterized by T p . Hence Weiss 
oscillations in magnetoconductivity in bilayer graphene are less damped with temperature 
compared to 2DEG system. 
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VII. CONCLUSIONS 



We have investigated the diffusive magnetoconductivity component a yy in bilayer 
graphene in the presence of a perpendicular magnetic field and a one-dimensional weak 
electric modulation. In this work, we focus on the Weiss oscillations in magnetoconductiv- 
ity. We have compared the results obtained with those of electrically modulated monolayer 
graphene as well as electrically modulated conventional 2DEG system. We find phase differ- 
ences of n and 2-7T in the magnetoconductivity oscillations compared to monolayer graphene 
and conventional 2DEG system which arises due to the chiral nature of quasiparticles in 
graphene. We also find that the oscillations due to modulation in the magnetoconductivity 
are enhanced and less damped with temperature compared to conventional 2DEG system 
whereas they are less robust with temperature compared to monolayer graphene. 
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